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To give a common theoretical description of liquid phases of the charged pion matter in a wide 
temperature interval, the relativistic quantum (p^ type model is considered. The liquid states of 
pion condensate and hot pion matter are investigated. 

PACS numbers: 21.30.Fe, 13.75.Lb, 03.75.Hh, 64.70.-p 



O 
O 

(N 

> 
O 



> 
o 

o 



o 



X 



INTRODUCTION 

The possibility to observe the pion condensate in 
heavy-ion collisions had served as a subject of theoret- 
ical and experimental investigations during last decades. 
This condensed state is often imaged like multipion 
droplets or/ and pion liquid. The pion condensate is also 
a necessary component of neutron stars. The theoretical 
description of pion condensate is usually given within the 
framework of phenomenological models leading to the cp'^ 
type self-interaction between pionic degrees of freedom 
(see, for example, [H, 0, 

Since the pions are the Goldstone bosons for sponta- 
neously broken chiral symmetry, a realistic description 
of pion subsystem can only be achieved on the basis of a 
model respecting chiral symmetry. However, if the pion 
matter is at the temperature much lower than the tem- 
perature of the chiral phase transion (T^ « 150 MeV), 
a chiral perturbation theory is applicable. Using per- 
turbation scheme, we would like to point out the gen- 
eral property of the chiral models, namely, in third- 
order approximation they result in an attractive two- 
body interaction (associated with ip'^ term) and three- 
body repulsive interaction (associated with (p^ term). To 
demonstrate it, we appeal to the Skyrme j^] and Wein- 
berg [5j models, where interactions can be presented as 
l^Skyrmc(;^) = -f^m^TviU + W- 2)/4 = f^m^iiyy2\ - 
^.4/4! + 1.6/6! - • • •) and lVci„bcrg(i') = fW^Vi"^ + 
= flm'^{v'^/2 - v^l% + iyV32 -...), respectively. 
Here we use notations for U — exp (zf • tt/Ztt) G SU{2), 

= if' I j"^; and m are pion decay constant and pion 
mass, respectively. Note that these interactions are the 
limited functions of v. Also note that the theoretical 
Gasser-Leutwyler interaction Q of order is attractive, 
and the experimental data justify the existence of three- 
pion correlations in heavy-ion collisions [t*!. Therefore, if 
we limit ourselves by consideration of the third-order ap- 
proximation (leading to the the relativistic quantum (p^ 
type model), we can expect to observe first-order phase 
transitions (PTs) in pion subsystem. 

Indeed, it is already known from molecular physics that 
the non-relativistic ip^ model, in contrast to the ip>^ one, 
allows us to observe not only second-order PT into con- 



densate state but also its gaseous and liquid phases. It 
is reached by means of solving the corresponding Gross- 
Pitaevskii equation [8,, J3J . Remark that the model has 
been already used in nuclear hydrodynamics and quan- 
tum field theory [lol | . Here we are trying to obtain similar 
results (without degrees of freedom of 7r° mesons) in the 
context of the physics of superdense ions and neutron 
stars, where pion condensate plays an imp_ortant role in 
the softening of nucleonic equation of state . Clearly, in 
this problem, we should take into account the essentially 
different nature of interactions in molecular physics and 
pion subsystem. Moreover, we would like to investigate 
the conditions of existence of liquid phase at high temper- 
atures (at the temperatures higher than the temperature 
Tcond of the second-order PT into condensate state) with 
the use of the same model. It is possible that the hot 
TT+TT^ liquid can be created in relativistic heavy-ion col- 
lisions, where the pions play a dominant role at the final 
stage of the reaction. Note that collective phenomena in 
particle-nucleus and nucleus-nucleus collisions are well- 
established and play significant role. An existence of the 
hot pion liquid has been already predicted in Ref. 



Since the liquid condensate of tt+tt" mesons thermo- 
dynamically differs from gaseous condensate, then such a 
difference should be taken into account, when attempts 
to register an appearance of a condensate in heavy-ion 
collisions are performed. On the other hand, the system 
of the large number of pions at high temperatures, as we 
shall see below, can also be in a liquid phase which can es- 
sentially affect on nucleonic dynamics. We say "liquid", 
when we deal with a dense phase in thermodynamics. 
However, the same "liquid" is a state with a high mag- 
nitude of pionic field from the field-theoretical point of 
view. As a consequence, amplification of pionic field can 
result in the creation of a proton-antiproton pair. 



THE PHASES OF PION CONDENSATE 



The Lagrangian density £(7r^, tt) of the model is 



£ = d^^Tx^d^-K - m^TX^-K + ^(^^7r)2 - |(^^7r)3, (1) 
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where the normal ordering of operator fields is assumed. 
We deal with the case, when electromagnetic interaction 
is neglected. Here it is also supposed that A — vn? /g'^, 
B = 3m^X/2g'^, m = 140 MeV, g and A are model pa- 
rameters which should be fitted. 

Constant g plays a role of pion decay constant rede- 
fined in medium. Model parameter A is introduced to ac- 
count effectively the higher-order terms of chiral interac- 
tion expansion. To analyze the range of A, let us present 
chiral interactions (considered above, for example) as 
V{v) = a2v'^ - aiv^ + a^v^^Xiv), where a„ = |y(")(0)|/n! 
and = 27rW/(7^. It turns out that \{v) is a smooth 
function, which is decreased from 1 to 0, when v runs 
from to oo. Thus, if even \u\ > 1, the constant A intro- 
duced instead of function should be less than 1 in 
order to relate our phenomcnological approach with the 
chiral theory. 

Our approach to this model is based on ansatz on appli- 
cability of the mean-field approximation (MFA) in a wide 
temperature interval. Moreover, we are limiting ourselves 
by two cases: i) T = 0, when the operator field tt can be 
replaced by a classical complex field (j), also called as the 
order parameter, describing a condensate; ii) T > Tcond, 
when there are no anomalous expectation values and the 
field TT coincides with quantum fiuctuations x- 
is to investigate the liquid-like states of 7r"'"7r~ matter in 
these two regimes. 

It is appropriate for us to begin from the investigation 
of different phases of the non-uniform pion condensate at 
T = 0, which appears in the neutron stars and the nuclei 
with density higher than saturation one. This situation 
is modeled by Lagrangian density In fact the 

replacement of tt by is analogous to the transition from 
quantum electrodynamics to the classical description of 
electromagnetism, when a big number of photons are in 
approximately the same state. In our case, the presence 
of a big number of pions in a single state (Bose-Einstein 
condensate) permits us to introduce the classical function 
(j). In a contrast to Maxwell theory, €{(/)*, (j)) contains the 
quantum constant h explicitly (?i = 1 in our units). On 
the other hand, the quantum meaning of classical fields 
is reviewed in Rcf. 

From the variational derivative of the corresponding 
classical action functional, one obtains the following evo- 
lution equation: 

^-V2^</. + m20-A|0|2^+|HV = O, (2) 

which serves as a relativistic generalization of the cor- 
responding Gross-Pitaevskii equation without external 
trap potential 

An applicability of Eq. ([2]) demands that the following 
conditions should be satisfied. First, the total number of 
pions (multiplicity) must be large enough because only 
in this case we are authorized to use the concept of Bose- 
Einstein condensation. Second, in order to replace the 



field operator by the classical field we have to assume 
both diluteness and the fact that the temperature is low 
enough. This allows us to ignore both the quantum and 
thermal depletion of the condensate. 

The solution of Eq. ^ essentially depends on physical 
parameters and boundary conditions. Here we shall deal 
with exact static quasi one-dimensional solution ("bub- 
ble" ) which can be presented as: 



0(t, z) — rik{z) exp (^-~imt\/ 1 — k^^ , 



(3) 



where rik{z) is a real function. 

The time dependence is in a phase, corresponding to 
the non- vanishing charge Q. Such a substitution results 
in a very large class of solutions rjkiz) in terms of ele- 
mentary and elliptic functions [is', 14]. However, we are 
interested in stable solutions arising from ordinary charge 
conservation. 

For the boundary conditions lim rik{z) = 0, 

z — ^±oo 

lim r]l{z) = 0, there is a solitary wave solution of the 

z—>±oo_Z 

form [i|]: 



2gk 



(4) 



Vl - 4Afc2cosh(2fcr 
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The result obtained describes a narrow layer of the dense 
phase trapped between dilute phases and leads to a lo- 
calized energy density. This model solution can describe 
a region of dense pion matter with respect to the axis of 
nucleus-nucleus collision. 

One can simply prove that the solution ((3]) is pseu- 
doscalar: it is necessary to replace the coordinate z and 
the (quasi)momentum k by —z and —k, respectively. As 
it must be, this operation changes a sign of (^(i, z). 

An order parameter in this system is the scalar density 
n(z) = \(l>{t, z)p. To determine the phase of pion matter 
in dense region, we introduce the charge Q (difference 
between the numbers of particles and antiparticles) and 
the total energy E per unit area: 



Q 

E 



8g ^ — arctanh-i 



Va 

Wl - k^Q 



1 - VI -4Afc2 
1 + V1-4AF' 



g^rak 



, (1 — 4Afc^)arctanh— 



VI - 4Afc2 
2fcVA 



(5) 



(6) 



These quantities are parameterized by k as independent 
variable. Also, it is not hard to see that, varying fc, the 
behavior of functions Q, E is mainly determined by pa- 
rameter A. At this time, the model parameter g influ- 
ences only the magnitude of these characteristics. 

Since the charge Q is conserved in time, there exists a 
non- vanishing chemical potential ^ = dE/dQ, which is 
calculated by the formula: fi = {dE / dk) / [dQ / dk). The 
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dependence of /i on can be explicitly found. However, 
we do not adduce it here because of its cumbersome form. 
Note only that > 0. It is in accordance with the result 
of 0], where multipion droplets are considered at low 
temperature. 

As shown in Fig. 1, the central density v = n{Q)/g^ 
and the chemical potential /i present backbendings typ- 
ical for the first-order PT. The transition point, given 
by the crossing point in e = E/mg^ versus q = Q/g^, 
corresponds to a Maxwell construction in the diagram 
of /X versus q. However, the system should never ex- 
plore the backbending part of the diagram because it is 
a metastable state. It is clear that dense phase is associ- 
ated with a liquid while dilute phase is a gas. We want 
to stress that both branches are quantum fluids. 

We see in Fig. 1 that the first-order PT in pion con- 
densate takes place at 0.25 < A < 0.268 and A < 1 as 
is argued above. Note that, at A « 0.268, the stable 
and metastable solutions coincide. It defines a critical 
point associated with a second-order PT: at this point the 
derivative of as a function of q diverges. For A = 0.25 
the attractive two-body interaction prevails and the sys- 
tem tends to collapse. In this case the maximum charge 
is limited by g w 5.3. Assuming an existence of the first- 
order phase transition, we put A = 0.26 in this paper. 
The value of constant g will be discussed in the case of 
hot pion matter. 



THE HOT PION LIQUID 

Now we focus on the description of the tt+tt" matter at 
high temperatures in MFA, when the field tt is described 
completely by quantum fluctuations x- We represent op- 
erator field X ill finite volume V as 



Up{x) 



exp {—lE-pt + ipx). 



(7) 



(8) 



The dispersion law is ~ y/v"^ + ml, where m, is ef- 
fective mass which should be found; aj^, b'^ (op, bp) are 
creation (annihilation) operators of particles and antipar- 
ticles, respectively. The commutation relations between 
operators are standard: [ap,ajj] = <^p,q, [&p-^q] = '^p,q- 

The Hamiltonian operator, H — J {'HQ + W)d'^x, based 
on C{x'^,x), is determined by the following terms: 



no = dtX^dtX + '^X^'^X + mlx^X, (9) 
W = {m'-ml)xh-^{xhr + %h?- (10) 
One has 



2! 



3! 



H = Y,Ep{np + np)+ J 



(11) 




FIG. 1: Central density v, total energy e, chemical potential 
jj, per mass m, in dimensionless units, as functions of the 
reduced charge q. 



where rip = a^ap and fip = b^bp are the number opera- 
tors of particles and antiparticles, respectively. 

Assuming that the interaction term is the perturba- 
tion, thermodynamic potential can be written as 



7.P 



_ „/3(7A'--Ep) 



VWmfa, 



(12) 



where fj, is the chemical potential, (3 — 1/T is inverse 
temperature, 7 = ±1. 

In our approximation, the interaction term is as follows 



Wmfa = {m^ - ml)a + m^{ba^ - 2aa^), 



(13) 



where a = (x^x) is the order parameter, a = l/2g'^, 
6== 3A/2g4. 

We find thcrmodynamically conjugate variables a and 
TO* from the extremizing conditions dQ/dm^ — 0, 
dQ/d(7 = 0, which result in the gap equation. 



and the mass ratio, 



1 



7,P 



2E^ e^(^p-w) - 1^ 



(14) 



ml/m^ = 1 - Aaa + 36ct^ 



Further investigation of thermodynamic properties of 
hot pion matter is carried out in the thermodynamic 



4 



limit, when the replacement of summation by integration 
over momentum takes place. Here we neglect a contribu- 
tion of the surface. 

It turns out that, to solve correctly the constraint, it 
is helpful to operate with dimensionless variable s in- 
stead of cr. It leads at once to the splitting of the 
value interval of a into two branches: a = {2g^/9X)[l ± 
y/l — (9/2)A(l — s^)]. Such a representation of cr gives 
us two equations F± (s, T, fi) — 0, where 



F±is,T,^i) = l±Jl-h{l-s 



A / 3m 
2 [2^ 



E 



dk 



1 



(15) 



The solution of these equations is the function s(T, /i) 
which is calculated numerically. 

One can say that this splitting reflects the existence 
of two phases of hot pion matter. If s = 1, we find 
that a single phase, corresponding to ideal gas, survives. 
Switching on an interaction (s 7^ 1), a new phase, called 
as hot pion liquid, can arise. The appearance of this 
phase is possible due to the sixth-order term. 
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FIG. 2: The dependence of the scalar density a (in MeV^) 
and the reduced effective mass s = m,/m on temperature. 

In order to compute thermodynamic functions, one 
needs to know the values of two parameters of our model, 
namely, g and A. We have already assumed that A = 0.26 
to achieve liquid-gas phase transition in condensate. To 
fit constant g, we try to appeal to experimental data 
obtained by DLS Collaboration in Berkeley and extract 
necessary information from dilepton spectra produced in 
proton- nucleus reaction 15| . If we believe that the main 
contribution to the spectrum comes from hot tt+tt" an- 
nihilation, one can conclude that the pion annihilation 
threshold is not equal to twice free pion mass (2rn) but 
smaller and approximately is 2ms = 260 MeV. Our fit 



to the spectra [l5|, which depends on the dielectron in- 
variant mass, gives us the temperature of the fireball, 
formed in this collision, equal to T = 75 MeV. There- 
fore, to find g at chemical equilibrium (^ = 0), we need 
to substitute these data in our equations. One obtains 
that 5 « 8 MeV from equation F+(s, T, 0) = 0. Using the 
equation F_(s,T, 0) = 0, another value of 17 is derived. 
However, the liquid-gas PT is not observed in this case. 

Assuming chemical equilibrium, there is the only pos- 
sibility to observe the transition into liquid phase with 
changing temperature. Indeed, the dense "liquid" phase 
appears with increasing temperature (see Fig. 2) that is 
in contrast to molecular physics, where the number of 
particles is conserved. This phenomenon has been al- 
ready pointed out in [llj. The critical temperature of 



the PT in our model is about 42.3 MeV and less than 



136 MeV as it was predicted in [ll|. However, we should 
conclude that, at T = 75 MeV of fireball, the pions are 
in the liquid phase, if the scenario, when 5 « 8 MeV, 
is realized in nature. This outcome demands additional 
theoretical and experimental verifications. 

Note that the calculations with non-vanishing chem- 
ical potential can be carried out and will be published 
elsewhere. 
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